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Abstract 

For Non-Relativistic QCD the velocity renormalization group correlates the renormalization 
scales for ultrasoft, potential and soft degrees of freedom. Here we discuss the renormalization 
of operators by ultrasoft gluons. We show that renormalization of soft vertices can induce new 
operators, and also present a procedure for correctly subtracting divergences in mixed potential- 
ultrasoft graphs. Our results affect the running of the spin-independent potentials in QCD. The 
change for the NNLL ti cross section near threshold is very small, being at the ^ 1% level and 
essentially independent of the energy. We also discuss implications for analyzing situations where 
mv^ ~ Aqcd- 
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I. INTRODUCTION 



The use of effective field theory techniques in non-relativistic fermion-antifermion systems 
has generated an encouraging record of successes. This framework has apphcations to de- 
scribing the dynamics of the bound quarkonium state as well as to heavy quark production 
and decay (for reviews see Ref. The formulation of a non-relativistic field theory for 
fermion-antifermion systems in QED and QCD has gone through several stages of develop- 
ment. In Refs. 0, a method was formulated for separating the short distance physics at 
the scale m from the long distance physics at the non-relativistic momentum and energy 
scales, mv and mv"^. Subsequent work served to clarify the power counting in u, and the de- 
scription of the low energy degrees of freedom. In particular the relevant low energy degrees 
of freedom have been classified into potential {E ~ mf^, p ~ mv), soft (i? ~ p ~ mv), and 
ultrasoft (E' ~ p ~ mv"^). 

It was realized in Ref. |^ that it is necessary to distinguish ultrasoft and soft contributions, 
and in Ref. that there was a problem with simultaneously power counting the ultrasoft 
and potential terms in the original NRQCD action. In fact, the ultrasoft gluons destroy 
the power counting in v unless their Lagrangian is multipole expanded @, ^j. This can 
be formulated in an elegant way by introducing more than one type of gluon field in the 
action [0]. Ref. p made the observation that since potential gluons do not propagate they 
can just as well be integrated out of the theory. In Ref. it was pointed out that in 
certain situations the matching coefficients can be efficiently computed from what is now 
known as the hard region of a diagram. This idea was formalized for other momentum 



regions and more general situations with the advent of the threshold expansion |[T^, which 



also emphasized the importance of soft momenta. Finally, the relevance of soft gluons for 



correctly running in the low energy theory was realized in Ref. |Tl |. 

In Refs. P, H^, [13| the authors went further to suggest that a string of effective theories 
could be defined by exploiting the hierarchy m ^ mv ^ mv"^. We will refer to the theory 
for scales mv < /i < m as mNRQCD , and for scales /x < mv as pNRQCD |^.^ A 
correct separation of the mv and mv^ scales is important since we would like to be confident 
about what effects can be reliably computed perturbatively, for instance in the case mv"^ ~ 
^QCD ^ ^v. The mNRQCD-pNRQCD setup appears to rely on being able to treat the 
energy and momentum scales of the non-relativistic quarks and their corresponding cutoff 
scales as independent. For static heavy quarks this is the case because the dynamics does 
not correlate the quark separation r ~ l/{mv) with energy fluctuations E -C mv. 

In Ref. [0 it was pointed out that for dynamic heavy quarks the non-relativistic disper- 



sion relation E = /{2m) couples the energy and momentum scales. The authors therefore 
proposed matching directly at the scale n = m onto a potential-like theory with both soft 
and ultrasoft degrees of freedom (referred to as vNRQCD). For scales fi < m the correlation 
of energy and momenta is accounted for since the ultrasoft and soft scales, fiu and /i5, are 
related, fiu = fJ^'^/m = mz/^. The running in the dimensionless velocity parameter from 
z/ = 1 to ^ of order the velocity of the two particle state, sums logs of both the mo- 
menta and the energy at the same time, and is referred to as the "velocity Renormalization 



Group" (vRGE) [|1^, |T5|. Within dimensional regularization the factors of /i^ or /i^ multi- 



plying operators in the renormalized effective Lagrangian are uniquely determined from the 



^ Note that we do not use the term NRQCD to refer to a theory for /i > mv. We reserve NRQCD as a 
generic label for the effective theory(s) which describe cc, bb, and ti bound state effects. 
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V power counting in d dimensions [0. This point is also discussed in Sect. |I|. 

With the vRGE, the running of potentials and currents was worked out in Refs. 



15], O, IT8| and applied to tt production near threshold in Refs. IH, |20|. The running 



of the static potential due to ultrasoft effects in the mNRQCD-pNRQCD formalism was 
computed in Ref. ||2T]]. In Ref. it was shown that the vRGE could be used to predict 
In a contributions in QED bound states, and for positronium the In^ a hyperfine splitting 
corrections and \v? a corrections to decay rates were reproduced, and the In^ a Lamb 
shift was predicted. In Ref. it was shown that the correlation of energy and momentum 
scales is necessary to compute QED corrections involving In'^ a with k > 2. More recently, 
the running of operators in the mNRQCD-pNRQCD framework were computed ||2^, and 
the original formalism was modified in Ref. to include the correlation of potential and 
ultrasoft cutoffs. 

The main purpose of this paper is to discuss two aspects of the ultrasoft renormalization of 
operators for mv'^ ^ Aqcd- We first point out a new set of soft operators which has zero tree- 
level matching, but is induced by mixing from ultrasoft renormalization of soft time-ordered 
products. These operators vanish for QED bound states such as Hydrogen or positronium. 
However they affect the leading logarithmic running of the two spin-independent 1 / QCD 
potentials as discussed in Sec. 0, and were not included in Ref. [|15|. We also compare our 



results to those in Ref. where analogous operators were included. In Ref. ||2J] two results 
were reported, one for scales mv < fi < m (mNRQCD ) and one for /i < mv (pNRQCD). For 
the potentials we disagree with the mNRQCD results because we find that for dynamic 
quarks there is no corresponding momentum region we can identify since the renormalization 
from ultrasoft gluons is always present. For the pNRQCD results we find agreement, however 
only if we force the mNRQCD-pNRQCD matching scale and the pNRQCD energy cutoffs 
to always be correlated. We also show that our approach reproduces the In^ a energy for 
muonic-Hydrogen with additional massless fermions p6[ . 

Second, we formulate a procedure for correctly subtracting ultrasoft divergences in di- 
agrams containing both ultrasoft and potential loops. Our results for the running of the 
l/(m|k|) and 1/k^ QCD potentials differ from Refs. |1T7| , |T^ , |2T1 , essentially because 
we find that there is an additional set of operators O^i, Od which must be renormalized. 
For matrix elements which do not involve additional divergences, the sum of these operators 
reduce to effective l/(m|k|) and 1/k^ potentials. For these effective potentials we agree with 



the n < mv pNRQCD results of Ref. |24] (again only if we demand that there is always a 



correlation between the mNRQCD-pNRQCD matching scale and energy cutoffs). However, 
our results also apply to matrix elements with additional divergences such as corrections to 
the QQ current correlators (7(0,0), where effective l/(m|k|) and 1/k^ potentials do not suf- 
fice. Numerically the change from the earlier QCD results is quite small, being < 1% for the 
normalization of the tt cross section. Our results have implications for the less perturbative 
situation mv"^ ~ Aqcd for dynamic quarks. They imply that NRQCD in this situation may 
be more non-perturbative than sometimes assumed. 

In Sec. |I| we review some of the formalism that we will need. We outline the issues 
associated with the ultrasoft renormalization of operators in Sect. |ITT|as well as our solutions. 
This is applied to muonic Hydrogen in Sec. 0, the QCD potentials in Sec. |V|, and the 

effect on the l/(m|k|) and 1/k^ QCD potentials is taken up in Sec.|VJ. We discuss the 
running of the production current in Sec. |Vli| , as well as implications for the scenario where 
mv"^ ~ Aqcd- In Sec. |V111| we discuss results for the production current correlators and give 
numerical results for the e~^e~ —>■ ti cross section near threshold. Conclusions are given in 
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Sec.m 



II. FORMALISM 



We begin be reviewing some aspects of the NRQCD formalism used here. As a way of 
motivation it is useful to recall that a properly constructed effective theory should satisfy 
the following requirements: 

1. reproduce the IR divergences of the full theory in its entire region of validity, 

2. have a well defined power counting (in v for our case), 

3. have no large logs in matching calculations, 

4. start with a regulator independent Lagrangian. 

The first property ensures that we have included the correct degrees of freedom. The need 
for the second property is obvious. The third property follows from the first, and the 
matching will be independent of the IR regulator as long as the same choice is made in 
the full and effective theories. The fourth property is necessary so that physical predictions 
are independent of the method used to regulate ultraviolet divergences. In practice we 
choose dimensional regularization which makes it easy to preserve the symmetries and power 
counting. For this case the fourth constraint dictates that the action is independent of d 
when expressed in terms of bare quantities. 

The physical system we wish to describe is that of a heavy fermion and antifermion with 
mass m, and energies E ~ mv^, and momenta p ~ mv in the cm. system where v <^1. The 
possible degrees of freedom include p|, H, ^, ^ |], ^ |], ||, [1^, |Tl| heavy potential quarks and 
antiquarks {ip^, Xp)i ultrasoft gluons, ghosts, and massless quarks (A^, c, v^^s), soft gluons, 
ghosts, and massless quarks (A^, Cg, ipq), potential gluons (A^ and soft heavy quarks and 
antiquarks (V'q, Xg)- Here, ultrasoft gluons are the gauge partners of momenta ~ mt;^, while 
soft gluons are the partners of momenta ~ mv. It is essential that we include both soft and 
ultrasoft gluons for all scales less than m. Evaluating a generic QCD scattering amplitude 
in the region of validity of the effective theory gives logarithms 

ln(E2), ln(p2), ln(k2), (1) 

where E denotes the cm. energy, p a quark momentum, and k the momentum transfer. Both 
ultrasoft and soft gluons are needed to reproduce the ln(£'^) and In(k^) terms. Furthermore, 
both types of logarithm are needed for all scales below m since both ultrasoft and soft 
running feed into anomalous dimensions induced by potential loops such as the production 
current and two-loop renormalization of 1/m? operators in QED p2|. Simultaneously 



including soft and ultrasoft gluons is also in agreement with the threshold expansion |T0[ 



where both soft and ultrasoft regions of energy and momenta are included in calculations in 
a separated form at scales ~ m. When these degrees of freedom are not treated separately, 
such as in the mNRQCD-pNRQCD approach for mv < fi < m, the v power counting breaks 
down, and a 1/m expansion with static quarks is used ^4\. However, it is argued that a 



power counting for dynamic quarks exists in pNRQCD for scales fi < mv |24] 
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To distinguish the mv and mv'^ scales we made use of a phase redefinition for the potential 
and soft fields 



Tl 



-ik-x , 



'k[.X) 



(2) 



where k denotes momenta ~ mv and d^(t)k{x) ~ mv (j)k{x)- Since here we will only be 
interested in cases with perturbative potentials (such as Aqcd ^ mv'^) we simplify the hst 
of degrees of freedom by integrating out potential gluon exchange and soft heavy fermions at 
the scale m following Ref. 
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This choice may not be unique, but does allow us to satisfy 
our effective theory criteria. For instance, in Ref. |jl6| it was shown that vNRQCD correctly 
reproduces all the infrared logs in QCD for four-quark Greens functions at one loop in its 
entire region of validity. 

The effective vNRQCD Lagrangian can be separated into ultrasoft, soft, and potential 
components, C = Cu + Cs + Cp. The presence of two types of gluons immediately brings up 
the issue of double counting. To avoid double counting the effective theory is constructed 
such that the ultrasoft gluons reproduce only the physical gluon poles where ~ k ~ mw^, 
while soft gluons give only those with ~ k ~ mv. The scales for the gluon momenta are 
influenced by the quark propagators, so the quark-gluon interactions must be constructed 
in such a way that we will not upset this scaling. In Cu this is achieved by the multipole 
expansion of interactions ^, which ensures that ultrasoft gluon momenta are always much 
smaller than the quark three- momenta. The ultrasoft Lagrangian is P, 

\2 



p 

+ ... 



IP 



iD) 



2m 



P 

8m^ 



x) \ - \g7GI 



(3) 



where G!^'^ is the ultrasoft fleld strength. In dimensional regularization the covariant deriva- 
tive has the form D'^ = +ifilj Qu^^, where = ttlv"^ and Qu = guifJ-u) is the renormalized 
ultrasoft QCD coupling. Note that the covariant derivative only contains the ultrasoft gluon 
fleld and that the ultrasoft Lagrangian has the form of the multipole-expanded HQET La- 
grangian. For convenience we suppress (throughout this paper) the renormalization Z factors 



that relate bare and renormalized quantities. The soft Lagrangian has terms |TJ, |T5 

1 
4 



16| 



(4) 



2 2e 



E 

p,p',q,q',cr 



+z^p,t [c,,, c,]Y^'^^ + (^p,t T^ZM V^p) (^,.7^T^^,) I + (^ ^ X, T ^ f ) , 

where G'^'^ is the soft gluon fleld strength and gs = fi's(/is) (/i5 = rrii') is the soft QCD 
coupling. The tensors uliu\ W^n) , z\j'\ and y^'^^ are functions of p', p, g, q' that are generated 
by integrating out soft heavy quarks and potential gluons. Their explicit form can be found 
in Ref. |17]. Finally the potential Lagrangian has terms il^, ITT 



r 

"n 



^(p, P') V^>pxLp,x-p + /i|^i^r'''(p, p')(^u/i^Ay) 



^l.T^V^pxLp,r^x-p 



+ 



(5) 
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where fi's^V and iig^F^^'~^ are functions involving bare Wilson coefficients. In the ffist term 
spin and color indices on V and the fermion fields have been suppressed. Matching pertur- 
batively at m and integrating out the potential gluons generates the terms 



r (p, p') = {T^ ® f^) 



(T) 



mlkl 2im?'k^ 



(T), 



k2 



+ 



+ 



V. 



(T) 



+ 



V 



(T) 



(T) 



A(p',p) + ^T(k) + ... 



+ (1® 1) 



V, 



(1) 



+ 



(Ti + (72 



2ivF^k, 



A(p',p) = 

i fABC 



. S-(p^xp) 

k2 



k^ ' m|k 
T(k) 



CTi - (72 



3 k • CTi k ■ cr2 



k4 



■r 



(6) 



where k = p' — p. The factors Vj{v) are renormalized Wilson coefficients, and the coefficient 



Vc in Fj"^^ is fixed by reparameterization invariance as shown in Ref . [0 . Some further 
operators are required to renormalize Cg^p and will be discussed later on. 

It is worth noting that the factors of /i^ and /i^ in Cu^s,p are uniquely determined by 
mass dimension and v power counting in ci = 4 — 2e dimensions ||16[. A scaling with v 
is assigned to the effective theory fields so that in the action their kinetic terms are order 



This gives ip^. 



[mv 



,3/2-. 



[mv 



2M- 



and Ai^ ~ 



Since ~ mu^, the 



[mv 

renormalized combination QuA^ must be multiplied by /z^ ~ {mvY so that this gluon term 
also scales as mv"^. For the potential fermion terms in Cs displayed in Eq. the soft and 
potential fields give a (mf )^~^'^ and the complete measure gives {mv)~'^~^'^'^ . The factor of 



'mvY''' is therefore required to cancel the v from the measure and fields. For the 



four-quark operator in Cp the quark fields give 



mv 



-4e 



and the complete measure gives 



[mv 



'^{mv)~ 



3+2e 



so a factor of /il" 



[mv 



\2e 



multiplying V{p, p') is required. The factors 



of and fXg for any other operator can also be determined in this way. 

The power counting of an arbitrary diagram is determined entirely by the powers of v 
assigned to its operators. Since the factors of //^ and fig are already determined, we can 
work in (i = 4 dimensions. In general a graph is order ~ v^ with [0 



5)VI'^ + {k ~ 4)V, 



(^) 



(7) 



Here Ng is the number of disconnected soft subgraphs if potential lines are erased, the vertex 
factors V^^'^ count the number of insertions of an operator scaling as f^. The factors V^^^ 
are for purely ultrasoft operators, V^^^^ are for operators which contain ipp or Xp but no soft 
fields, and V^^'''' are for operators with at least one soft field. In this paper we will refer to 
the order in v of diagrams and operators as their value of 

6' = 6-5, (8) 



since this quantity is essentially additive with multiple insertions of operators [|22|. For 
instance, using the value of 6' for the power counting a T-product of the potentials Vi ~ f " 
and V2 ~ v^ scales as T[ViV2] ~ v"''^''. With this counting the Coulomb potential scales as 
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(a) 




FIG. 1: Examples of mixed soft-ultrasoft graphs with ultrasoft gluons which could renormalize 
Cs in Feynman gauge. The zigzag lines denote soft gluons, quarks, or ghosts. 

v~^, the l/|k| potential is order f° and the potentials are order v^. For later use we 

also define 

a = k-4: (9) 

as the order in v for which a soft operator contributes to the power counting formula in 
Eq. (0). 



III. ULTRASOFT RENORMALIZATION OF OPERATORS 

Since the effective theory simultaneously involves soft and ultrasoft gluons the question 
naturally arises in what manner the two types of gluons renormalize operators. As far as 
renormalization of is concerned the ultrasoft renormalization is equivalent to that in 
HQET. This is because after the equation of motion E = p^/(2m) is applied, the quark 
propagators in loop graphs in the single fermion sector become static 



£; + fco - p2/(2m) +ie P + ie 



(10) 



Here {E, p) are the energy and momentum of the external quark line. Thus, the renormalized 
Cu is exactly of the form of the multipole-expanded renormalized HQET Lagrangian. 

The ultrasoft renormalization of interactions generated by Cg and Cp is more subtle. The 
basic complication is that one must account for the fact that ultrasoft gluons renormalize 
operators which are non-local with respect to the mv scale (but are local relative to mv"^). 

First consider the renormalization of interactions generated by Directly renormal- 
izing Cs by ultrasoft gluons leads to gauge dependent Wilson coefficients as observed in 
Ref. To see this consider the graphs in Fig. |T] which involve ultrasoft gluons. In 

Coulomb gauge these graphs are identically zero. In Feynman gauge the sum of Fig. ||a 
and the quark wavefunction renormalization give a Cyi/e term which could contribute to 

the renormalization of U^J in Cs- Non-Abelian graphs such as the one in Fig. |l|b do not 
cancel this contribution. For Fig. |l|b, the intermediate soft gluon propagator has momentum 
(g + kY = q'^ + 2k ■ q + k"^ , but we cannot set the offshellness ~ {mv)"^ of the soft gluon to 
zero without violating the power counting and risking a double counting. Keeping ^ 0, 
Fig. 0b evaluates to zero and it is not surprising that the total result is gauge dependent 
since part of the calculation was done offshell. This issue can be avoided by only considering 
the ultrasoft renormalization of operators which can contribute as soft color sing lets 



such as time-ordered products involving two or more soft vertices and quarks and antiquarks, 
T{Cf C^ . . .). These products appear local as far as the ultrasoft gluon is concerned, and it 
is only these products which affect observables. Having ultrasoft renormalization only for 
these operators also avoids the predicament of having both ultrasoft and soft gluons in the 
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(a) 




FIG. 2: The zigzag lines denote soft gluons, quarks, or ghosts. Graphs (a) and (b) are examples 
of mixed soft-ultrasoft graphs, while (c) denotes an operator for soft Compton scattering off a 
potential. 

single heavy quark sector. Thus, we do not consider ultrasoft renormalization of a single Cs 
term. 

However this does not imply that the first ultrasoft renormalization of a graph involving 
soft vertices are for two loops graphs like the one in Fig. The reason is that we must 
first consider the case where only one pair of soft gluon fields in the time-ordered product 
are contracted. So at one-loop we must consider the renormalization of graphs such as the 
one shown in Fig. 0b. The ultraviolet divergences in these graphs need to be canceled by 
counterterms for the 6-field operators, shown in Fig. §c. For quarks, gluons, and ghosts 
these operators have the structure 

o'^i = 9t4' i^l' rS:i ^p) (xlp. x-p) (^-, rj^) v,) , 

= gt ri^i ^p) (xLp, v'aI x-p) {At, vf^^^ a^,) , (ii) 
o^^} = gtf^t v^p) rg x-p) (c-.r^c,), 

where color indices are suppressed and the factor of is determined by the procedure in 
Sect. |T[ Here the Fj are matrices in spin and/or color space and can be functions of p', p, 
and q^. This soft momentum dependence is identical to that in the graph on the LHS of 
Fig. ^ once the equations of motion are applied (e.g. = 0, ^(pg = 0). The superscript 
a > has been used in Eq. (P) and denotes the order in v for which these soft operators 
contribute to the power counting in Eq. (^. For our purposes the a = 0, 2 operators will be 
sufficient. The operators in Eq. ([TT|) also have Wilson coefficients C2°^''(z/). The tricky thing 
about the operators is that the tree level matching onto their Wilson coefficients is 
zero, C2'^\i' = 1) = 0. An example of this matching calculation is shown in Fig. |^ for . 
At u = 1 the full theory graph on the left is exactly canceled by the time ordered product of 
soft vertices. Therefore, we get zero for the Wilson coefficient of the operator on the right. 
However, the ultrasoft loop graph in Fig. §b gives O^^^ a non-zero anomalous dimension 

(2) 

SO that for v < 1 the coefficient evolves and (z^ < 1) 7^ 0. The same counterterm is 
necessary to renormalize the divergences in Fig Furthermore, the operators O^i affect 
the anomalous dimensions for the 1/m? coefficients V2 and Vr in Eq. (j^), and were not 
included in Ref. This occurs through an ultraviolet divergence in the graph in Fig. ^ 
which needs to be canceled by V2,r counterterms. In Sees. |rV| and |V| we will give two examples 
of the implications of the operators in Eq. (|TTp. 

Let us now discuss the ultrasoft renormalization for two-body interactions generated by 
Cp. In general in graphs with two heavy quarks the equations of motion do not make the 
propagators static. For example, in the order two-loop graph in Fig. ^ the ultrasoft loop 
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FIG. 3: Example of the matching calculation for 021^ - Here the zig-zag lines denote soft massless 
fermions. At the high scale {v = 1) the graphs on the left exactly cancel so the coefficient of the 
operator on the right is zero. 




FIG. 4: Example of renormalization of two-body potential graphs by an ultrasoft gluon. The 
ultrasoft couplings are the order term from Eq. @ and the order term p • A in Eq. (|3|). 



momentum £^ is routed through an internal fermion line giving a propagator 

i = i (12) 

E + £0 + gO - qV(2m) + + gO - (q2-p2)/(2m) + ' ^ ^ 

where q = {q^, q) is the potential loop momentum. The multipole expansion has resulted in 
factors of the three-momentum £ being dropped, however the remaining propagator is not 
static in general since q^ 7^ p^. In fact, the ultrasoft loop induces an UV divergence after the 
dq^ and d'^i integrals have been performed, but before the sum over indices q is carried out. 



In Refs. [0, |T^] the subtraction of the ultrasoft UV divergences was made after carrying 
out the sum over q. Instead, the corresponding operators used to subtract the divergences 
should involve a sum over q. For example, for the sum of all scattering diagrams at order 
a^f " with one ultrasoft gluon the operators needed to subtract the ultrasoft UV divergences 
have the form 



^k2 



E (A + f-) [<'T'i^.xl,^T\.,] , (13) 



m 

p,p',q 



where the are functions of p, p' and q that will be given in Sec. ^I[ These operators are 
denoted graphically by the diagram in Fig. ^d, since they are essentially like the product 
of two potentials summed over the intermediate 3-momentum q. A similar set of operators 
exists for three-loop diagrams with an ultrasoft gluon at order To evaluate matrix 

elements of these operators using dimensional regularization we combine the sums with the 
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integration over the ultrasoft spacetime coordinate x 



, which leads to 



(14) 



Unhke the sums over labels on the fields, the free sum over q scales as v'"^"^, and contributes 
in determining the factor of /x^ in Eq. (|TB|). 

At the level of Fig. ^ it still appears ambiguous whether the sum over q in Eq. (|13D 
needs to be carried out. The key point is that the original ultrasoft divergent loop in Fig. |a 
acts like a one-loop sub divergence, despite the fact that it shows up at the level of the two 
loop graph. Therefore, it renormalizes each term in the sum of the operators in Eq. (|13|), 
rather than a potential V^(p, p') with the sum over q already carried out. 

This renormalization prescription can be illustrated by considering a graph with an ad- 
ditional ultrasoft fermion bubble as in Fig. To renormalize Fig. ^ we require both the 
1/e counterterm for the fermion bubble as in Fig. and a 1/e^ counterterm from Fig. 
which has an identical momentum structure to the operators in Fig. ^d. In dimensional reg- 
ularization the complete result for the sum of the three graphs has the divergent structure 



(15) 



where in the limit e — > all possible sub divergences ln(i?^//i^)/e and ln(k^//i|)/e are 
canceled by the counterterm contributions. If we only needed counterterms to be added 
after the sum on q was carried out then the sum of Figs. He and 0d 



/4 

k2 



(16) 



should only contain an overall analytic divergence. However, Eq. (^) has a non-analytic 
divergence. We note that as for the operators the Wilson coefficients Vki{i^) of the 
operators O^i vanish at the hard scale, i.e. Vki{i^ = 1) = 0. However, divergences from 
graphs such as those in Fig. ^ lead to a non-zero anomalous dimension so that Vki{i^ < 1) 7^ 0. 
Furthermore we emphasize that the renormalization of these operators is not equivalent to 
the renormalization of the QCD l/|k| and l/k^ potentials or to the renormalization of 1/r 
or l/r^ potentials. This will be discussed further in Sec. pl[ 



IV. MUONIC-HYDROGEN WITH MASSLESS ELECTRONS 

Our first example is a simplified toy model for muonic-Hydrogen with proton mass rrip — *• 
00, muon mass m fixed, and nj massless "electrons". In this theory the running of the 
coupling is simplified because there are no gauge boson self interactions. Furthermore at 
the order we are working the soft graphs simply involve massless leptons. This model was 
proposed in Refs. p4|, ^ as a test of the anomalous dimensions of the l/m^ operators (see 
also Ref. ||26| and reference one in H]). In this section we show that vNRQCD correctly 
reproduces the n/a^ln^ a energy levels in this toy model and is therefore in agreement with 
the QED limit of Pachuki's results in Ref. ||28|| . 

For this toy model we have {T^ (8) T^) (—1), (1 ® 1) (+1) ^^^d therefore define 

(1) (T) 

Vi = ^ -Vi- For the soft spin-independent vertices with rrip = 00 we only require 
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4°^ = l/k2 and Z^^^ = -CD{mi^)/i8m^) from Ref. |T5|. Here coifi) is the QED limit of the 



gauge invariant Darwin coefficient computed in Ref. [p9[ 



1 16 , 
1 + ^In 

3/5o 



am 



(17) 



where /3o = — 4n//3. The operator O^^ has the form 



0(0) 



(18) 



where 



r?Hg,P,p') 



(2^070 + k ■ 7) ^ (2^070 - k ■ 7) 



k2 + 2k ■ q 

r(o) = _] 



k2 - 2k ■ q 



(19) 



and k = p' — p. From tree level matching we have C2^J{1) = 0. Furthermore, since ultrasoft 

photons bring an extra f^, there is no anomalous dimension induced for this coefficient, so 

in fact C^J{i') = and this operator does not need to be considered further. The operator 
(2) 

has the form 



(2) 
2ip 



a 2 2w 



(20) 



From tree level matching one finds that its coefficient vanishes also at the hard scale, 
Cg)(l) = 0, but it gets a non-trivial anomalous dimension from the UV divergences in 
the ultrasoft graph in Fig. ^ plus wavefunction counterterms. 



„ 1 k^ a(mz/2) ,.„(o)\ 
+ w.fn. = ^ {tO^^J) 



2^ /4e27 



1 a;(mz/2) 



{tO?^),e2,, (21) 



where the soft vertices symbolized by the dots stand for the coupling Together with 

the respective soft divergence induced by the pull-up mechanism [jl8| this gives 



1 

TT 



(^^2t^)4e27 



(22) 



We see that the result would vanish if the two couplings were evaluated at the same scale. 
The result in Eq. (^) is canceled by a counterterm for C2J which, using the vRGE, gives 
the anomalous dimension 



77— r^') 



TT 



2a(mz/ ) — a(miy) 



(23) 



There is another possible contribution to the anomalous dimension in Eq. (|23|) coming from 



(2) 

the soft coupling renormalization graphs for the operator 02^, however these are exactly 
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FIG. 5: a) ultrasoft gluon graphs with p • A vertices (with wavefunction renormahzation on the 



other hne understood), b) soft graph with Zg'^ and Z^" ^' vertices, c) soft graph involving 02^5 . 
The zig-zag lines here denote massless soft fermions. 



{a'=2) 



l(2) 



canceled by the charge counterterm contribution associated to the factor in Eq. 
Similarly, an UV divergence in the analogue of Fig. 0b having only soft gluons is exactly 
canceled by a counterterm associated to the Zq term in the soft Lagrangian. Solving 
Eq. ( p3|) with the boundary condition C2\^{1) = gives 



4 

Jo 



In 



a;(mz/^) 


3 - 




~ 4 . 



CD[mi' ) — CDKmu ) 



(24) 



The anomalous dimension for V2 can now be computed from the graphs in Fig. ^ which 
are order v^e^. The graphs in Fig. ^d, |^c give the following three terms 



v-^V2 = — - aimv) a{m'u'^) + — c/)(mz/) [a(mz/)]^ H — ^ 
Of 3 3 9 



C^;^(z.)[a(mz.)]^ (25) 



Using the solution in Eq. ( p4|) this becomes 



4v, 



— a(mz/) a{mv'^) H — - CDimu'^) [a(mz/)]^ . 
3 3 



Solving this equation with the boundary condition V2(l) = iTa{m)/2 we find 

IT 



Cni'mi' ) a(mz/) 



The result differs from the coefficient of the four quark operator in Ref 
{mv < fi < m). Our result is in agreement with the pNRQCD result in Ref 



(26) 
(27) 

in mNRQCD 



TT 



air 



(28) 



but only if we correlate the mNRQCD-pNRQCD matching scale 1/r with the pNRQCD 
renormahzation scale Uus for the energy by enforcing Uus = mv"^ and 1/r = mv. We also 
agree with Pachucki for the ri/a^ln^a energy levels for this toy model. In Ref. |2^ 
a different expression for V2(z^) was inferred for the vNRQCD approach, however this is 
because the contribution from the graph in Fig. ^ was missing. 



V. QCD RESULTS FOR l/m^ POTENTIALS 



Our second example of the effect of the operators in Eq. (|TT]) is non-relativistic QCD 
for equal mass heavy fermions. In this case the soft degrees of freedom include soft gluons, 
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quarks, and ghosts. 



The operators , 

r(0) 



(0) Qm 



2A 



and of} 



in terms of the soft vertices U]nJ , Y'^^\ and Z^^^ from Ref. 



can be written in a compact form 

]], as summarized in Appendix 

A. Closing the respective two soft (gluon, quark and ghost) hues one obtains for the one- 
loop four quark matrix element a structure that is identical to the one-loop time-ordered 
product of the soft vertices ufj, Y'^^\ and Z^^^ [up to non-trivial Wilson coefficients which 
are suppressed in this equality]: 




—I 



T^®T^ + ... 



(29) 



0=0 



Just as in the matching calculation for the QED toy model the momentum structure of 
the soft time-ordered products agrees with the full theory result, so at z/ = 1 we find 
— 0. Similarly, since ultrasoft gluons bring an extra factor there is again no 



anomalous dimension induced for the operators and we have 



identically. 



For 0" = 2 it is necessary to keep the different color structures between the heavy fermions, 

(2) (T) 

SO in general we have to consider two types of operators, O^i with the fermionic struc- 
ture ['?/'p,T^-?/^pxLp/^'^X-p] and Of^''^^'^ with the fermionic structure [V'p/'^/'pX-p'X-p]- These 

operators are defined by Eq. ( p2D with Og^'*'^^'"^^ defined in Appendix A. For our purposes it 
is sufficient to consider only operators where the soft lines are closed in color space, and the 
operators Ogi'*^^'^^ Of^^^''^^ are defined with this convention. This is sufficient because 
for the renormalization of the potentials all external soft lines are contracted. From tree 
level matching one obtains that the coefficients of these operators vanish at the hard scale 



Cf^''^^\l) = 0, but as in the QED case they have a non-vanishing anomalous 



c£''"-'(i) 

dimension due to the exchange of ultrasoft gluons in diagrams of order a'^v^ as shown in 
Fig. Now since we are dealing with equal mass fermions, the p ■ A ultrasoft gluons can 
attach to any of the heavy fermion lines. Including the permutations of possible ultrasoft 
attachments we find 



+ 



+ . . . 



3 

Ca 



AQ2i 



V 3 12 4 



3m^ 



(30) 



where the ellipses denote other attachments including wavefunction diagrams. For SU(A^c 
QCD the color coefficients that appear here are 



- I 1 

Ca = Nc, Cf = — , Cd = 8Ci? — 3Ca, Ci = -CfCa 



(31) 
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The result in Eq. (|30D renormalizes the operators 



(2),(1) 
2a 



— \^ 



(0),(1) 



o. 



(2),(T) 
2b 



(0),(T) 
2i 



A^(2),(T) 
^2c 



(0),(T) 
2i 



(32) 



whose coefficients C*-^'' 



2a 5 



C'^^ , and Cjc'' vanish at the matching scale v = 1. Adding to 
Eq. (^) the soft divergences induced by the pull-up mechanism one arrives at the following 
anomalous dimension for the coefficients: 



.(2) 



(2) 



d 

lu2{j^) 



-2Ci 



u—C. 



(2) 



These equations are quite similar to Eq. 



Gtt 



7«2(z^) 



ICa 



iu2{y) 



(33) 



3D and their solutions are 



3(3o 



ln(w) 



CSV) 



3CA~Cd—4:Cf 
Wo 



ln(w) 



where 



w 



(34) 



(35) 



Note that as in the QED case there are additional possible contributions to the anomalous 

(2^ (\ T) 

dimensions of the operators ' which are, however, exactly canceled by pre-existing 

counterterms. This includes for example the purely soft diagram analogous to Fig. and 
the soft coupling renormalization graphs. The divergences of the soft diagrams are canceled 
by counterterms for U^J , Y^'^^ and the divergences of the soft coupling graphs by the 
counterterm associated to the factor g'^ in Eq. (|lT]). 

Finally, we consider the anomalous dimensions for the spin-independent potential coeffi- 
cients V2^f"\jy)- The graphs required for this computation are shown in Fig. The graphs 
in Fig. ^ and IpD were computed in Ref. |15[, however the contributions from the graph 
in Fig. ^ which are proportional to the coefficients were not included. Adding the 
contributions from all three graphs the results are 



z/— V(^) 
du ' 



32 



2(/3o + 7:Ca) oil{mv) + "^Ca as{mv)as{mi''^) — 4/3o Cfcii^) Cis{^y) 



d 



(T) 



/3o 
2 



1 + CdM 



24(z/) 



6 



28 - IIcdM + 26cfM 



+ ^{4Cf + Cd- 3Ca) as{mv)as{mv') - 2/3oCgV) as{muf , 



7C, 



14 



16 



Ci ag(mz/) — 5~ Ci <ys{mi')as{mi''^) —2(3o C^J{i^) asimuY . 



(36) 
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FIG. 6: a) ultrasoft gluon graphs with p • A vertices on all permutations of the four fermion lines, 
b) soft graph with Zq'^''^ uj^u''^ ■* or vertices such that a -\- a' = 2, c) soft graph involving 

one of the O^^ oj^^ O^J operators. 



The gauge invariant HQET coefficients that appear here are 

-20 32Cf 



-^Ca//3o 



^-2Ca/Po _j_ 



61L.F \ V. _ -13Ca/(6/3o)1 



.13 13C^. 

where z = as{mu) / asijn). The solutions for the coefficients in Eq. (^) are 

327rC^ 



47r aJm) z 



3/?o 



as{m) z ln(w) 



(37) 



(38) 



+ 



TT [C^(48Cjr + UCd + AC a) - /3o(33Ca + 32C^)] 



as[m) [z - Ij 



13C^(6/?o - 13Ca) 



+ a.(m) [z-^^-/^" - 1] - ^^^^^^^^^^^«.(m)zln(^) 



(/5o - 2C^ 
47rCi 



3/3o 



/ \ /-, \ ISttCi / \ 1 / n 
Q;s(m) (1 — z) -\ — — as{m)zm[w). 



For the color singlet channel V^'^^ = V^^^ — CpV^'^^ and the above results give 

8Ca 



Vi"\u) = TTCpasim) (z-l) 



-At: Cp «s("^) -2 



3/3o 



ln(w) 



(39) 



33 32^^ 90 A lOOC^ 



13 13^^ 13/^0 13/?o 
SttCf (3/?o - 1 IC^a) (5Ca + SCf) 



13C^(6/?o-13Ca) 



7rC^(/jo-5CA) 
{(^o-2Ca) 



■ajm z 



,1-2Ca//3o 



1] 



16nCF{CA-2CF) 
Wo 



as{m) zln(w) . 



The results in Eqs. (^) and (|39D are valid for all values of u between 1 and a vq of order 
the physical velocity of the quarks. This accounts for energy scales between m and mv'^ 
and momentum scales between m and mv, so our results apply to this entire region. Our 
results do not agree with the running of the 4-quark operators dgs, dgv, dvsi d^v obtained in 
Ref. |2^ for mv < n < m (mNRQCD ). This is because in Ref. ||2^ ultrasoft gluons are not 
included in the results for > mv. To translate the results in Ref. |2^ for scales < mv we 
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FIG. 7: Ultrasoft graphs for the renormahzation of the operators Of,-^ (left-right and up-down 
symmetric graphs are implied). 



note that V, 



is) 



-A-kCf D\ 



(2) 



TrCpiDfl — D)^l), and we again enforced a correlation 



l(2) 



by demanding that the pNRQCD renormahzation scale Vus = mp^, and the mNRQCD- 



pNRQCD matching scale 1/r = mv. With these restrictions our results in Eq. ( p9D agree 
with Ref. |2|] for the case ^ < mv. 

Note that there are also constraints on the relation between the cutoff for the momentum 
p of the quarks and the cutoff for the momentum transfers k = p' — p. Reproducing the 
known a^ln^a hyperfine splitting for positronium [3D|, requires that the cutoffs for these 
scales are correlated, since the calculation depends on simultaneously integrating anomalous 
dimensions that arise from soft and potential loops 

The remaining coefficients of the 1 / potentials are spin-dependent and not affected by 
the operator in Fig. ^c. They were computed in Ref. |jl5|: 



2tx 



(2C^-/?o) 
— —aAm] z 



a.Am 



(40) 



{\-2CaIM 



These expressions for and agree with Ref. |^ (however Vs disagrees) 
results for the coefficients in Eq. (HOf) were confirmed in Ref. 



The 



VI. QCD RESULTS FOR OPERATORS 0^, Od AND l/|k|, l/k^ POTENTIALS 

At order a^v^ there are three diagrams which have ultrasoft gluons (shown in Figs. ^,b,c) 
whose divergences are not canceled by counterterms from V2,r Jl^l- pointed out in 
Sec. |T| these diagrams should be renormalized by operators that do not have the sum over 
intermediate potential quark 3-momenta q carried out, as pictured in Fig. |^. Evaluating 
Figs. 0a,b,c we find that the required operators have the form 



p,p':q 

= -^^^^^ E + f^) H'T^^r>xl,,T\,,] , (41) 
p,p ,q 



(1) (1) (T) (T) 

giving the contribution ACp = Vki O k2 ~^^k2 ^ ki vNRQCD Lagrangian with Wilson 

(1) (T) 

coefficients V^/ and , respectively. The fig factors in Eq. ( ^T] ) are determined as in 
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sections ||,|T3. In Eq. ( ^T]) the functions /« are 



/o 

/2 



p • iq - pj 



(q_p)4(q_p/)2 



IP ^ P 



fi 



q ■ iq - pj 



(q - P)^ (q - P')' 



+ (P ^ P' 



P' 



iq-pj 



(q_p)4(q_p/)4 



(q^ 



P'V2 - p72) . 



(42) 



If the finite sums over q were carried out as integrals d^q or d"'q with n = d — 1 ^ 3 then 
the operators O^/ would reduce to l/|k| potentials. However, in general it is the operators 
in Eq. (^iD which are the fundamental quantities as discussed in Sec. |T|. The counterterms 
required to cancel the divergences in the two-body ultrasoft graphs in Fig. |^a,b,c plus the 
corresponding soft divergences determined by the pull-up mechanism [|l^ are 



2CaCi 
3^ 



(T) _ CAjCA + Cd) laJjHj) asifJ's) 



k2 



• (43) 



Next consider the order a'^f ~^ ultrasoft graphs. At this order there are two diagrams 
(shown in Figs. ^,b) with ultrasoft gluons which contain divergences that are not canceled 
by counterterms from V2,r,fei,fc2 H^- ^ P^^^ ^^e UV divergences in Fig. §3 is canceled by 
a ^^ki,k2 counterterm graph involving the function /2 as shown in Fig. ||c. Similar to the 
graphs in Fig. |^ the remaining divergences are subtracted by operators with sums over q, q' 
as pictured in Fig. |^. The required operators are 

p,p',q,q' 
p>p'>q.q' 



'c2 



(44) 



p.psq.q 



(1) (T) (T) 

where V^i , V^2 ■> '^cz corresponding Wilson coefficients that appear in the resulting 



V]:\>0).\> + V).^>0'^> + V^],^Og\ The /x^ 



contribution to the vNRQCD Lagrangian, A£p = V^i^Cji^ + VjpOiP 
factors in Eq. ( PD are determined as discussed in sections pi||rn[ The functions /io,i have 
the form 



ho 
hi 



[q - p j ■ (q - pj 



(q-pj^q 



q')'(q' 



p 



/^4 



iq - q j ■ iq - pj 



(q_p)4(q_q/)4(q/ 



P 



'^2 



+ (p ^ p ,q ^ q 



The counterterms needed to cancel the ultrasoft divergences in Figs, 
divergences determined by the pull-up [|l8l are 



(45) 

,b,c plus the soft 



5V, 



(1) 

cl 



CaCi{Ca + Cd) \as{^-u) oisi^is) 



VItx 



5V. 



(T) _ ^CaCx [as(/^, 



c2 



'c3 



2Ca 



e 



e e 

32 
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(46) 




FIG. 8: Graphs for the renormahzation of the operator which leads to the potentiah 



From the counterterms in Eqs. (43) and (46) the vRGE gives the anomalous dimensions 



-7«2(z^) , 



d_ 
du 
d 



(1) 



37r 

CaCi{Ca + Cd 



(T) 
k2 



(T) 



4CU 
37r 



Gtt 



-7«2(«^) 



3^ 
SCaCi 



-7«2(z^) 



37r 



7«2('^) 



{Ca + Cg 

32 



7«2('^) • 



(47) 



Note that possible soft graphs which could contribute to these anomalous dimensions are 
exactly canceled by counterterms similar to the QED example in Sec.|IT|. Using the 



boundary conditions V^^'"* (1) 



V, 



fci ^ 



ki 



ln(w) 



we find the solutions 

2Ca{Ca 



v«(.) 



3/3o 

CaCi{Ca + Cd) 



Cd 



V, 



c3 



'Wo 



3Po 
Ci + 



V 



c2 



3/9o 
16CaCi 

3/3o 



ln(w) 



ln(w) 



{CA + Cdf 

32 



ln(w) 



(48) 



where w is given in Eq. (|35|) . These are our final results for the Wilson coefficients of the 
new operators in Eqs. (^If ) and (|4^). The corresponding color singlet channel coefficients are 
given by V''*-' = V*-^^ — CpV^^^ (where coefficients for color structures not shown in Eqs. (^T]) 
and (|4^) such as are zero at this order). 

With the operators in Eqs. (|4T|) and (^41) all divergences in the usoft diagrams in 
Figs. |3a,b,c and |a,b,c are canceled completely. Therefore, the corresponding contribu- 
tions in the anomalous dimensions for the potential coefficients ^ in Eq. (^ obtained 
in Refs. [|T7i |l^ should be removed. The only remaining divergences that must be canceled 

(1 T) 

by ^ are from purely soft diagrams and are associated with the known running of the 
strong coupling a^. Thus we find 



87r 



Cd) [as{mv)f 



4v-i^) 



atimi') 



(An) 



where f3o,i,2 are the coefficients of the QCD beta function (in the MS scheme for ^2) 
the matching conditions at z/ = 1 [ff6l these equations give the solutions 



(49) 
With 



V, 



i7CA-Cd 



a^imu), v;^\u) = ^a^{mu), vP{u) = Arcafimu) . (50) 



Ci 
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Here oii\mv) is the QCD coupling with 3-loop running. 

It must be noted that in general the operators Okj and Ccj, are not directly related 
to the potentials l/|k| and in Eq. (^. The reason is that the sums over q, q' must 
be regulated in the same way as sums over p, p'. This has important implications for the 
cancelation of sub divergences as discussed in Sec. ^ and therefore also affects renormalized 
matrix elements as we will discuss further in the next two sections. However, if we take 
finite matrix elements of the l/|k| and 1/k^ potentials and operators Okj and Ocj-, then we 
can use E^/o = Eq/i/3 = /2 = l/(16|k|) and ^^o = Y.^,^>h^/2 = l/iGAn^k^). 

The result for the finite matrix element is then equivalent to the matrix element of effective 
l/|k| and 1/k^ potentials with modified Wilson coefficients: 



i% = V«(z.) + 6«^(mz.)V«(z.), Vfi = Vf)(z.)+4a,2(mz.)v(?( 



V.% = VPi'^) , VS = VPi'^) + 7ra',{mu) [vj^(z.) + 2V,T(^)] • (51) 

Substituting in the solutions in Eqs. and ( |5I1D this gives the following effective color 
singlet coefficients 

n,es = -^(Cf-2Ca) a,(mz/) + — a,(mz/) \n{w) , 

V^J^ = -AnC^af^ (mu) + ^^^^ aHmu) Hw) . (52) 

We emphasize that these effective coefficients can only be used for finite matrix elements, and 
in this case they should be viewed simply as a shorthand way of accounting for insertions of 
the Oki and Od operators and the l/|k| and 1/k^ potentials. The results for these effective 
coefficients agree with the results for the coefficients of the l/|k| and 1/k^ potentials in 



Refs. ||2l| , for scales < mv (pNRQCD), but only if we impose the scale correlation 
and 1/r = mu as mentioned previously in Sec.0. Again, our results do not 



agree with Refs. [21, 2^] for scales mv < jj, < m (mNRQCD). Furthermore, except for finite 
matrix elements, our results disagree with Refs. [|1^, [T8|, ^ since we have found that the 
anomalous dimensions are associated to a different set of operators which includes Oki-, Ok2, 
Oci, Oc2, and Ocs rather than just the l/|k| and 1/k^ potentials. In general these operators 
contribute in a different way for matrix elements with divergences. 

VII. RESULTS FOR THE PRODUCTION CURRENT 

The leading order production currents in the effective theory are of order and produce 
a qq pair in a ^Si or state: 

Ji,P = ^i'^(^^2)x1p, Jo,p = ^l{i(T2)X-p- (53) 

where Ci(z/) and Co(i^) are the corresponding Wilson coefficients. The fields are evaluated at 
spacetime coordinate x. The two-loop matching for Ci(l) and Co(l) were first considered in 



Refs. ||3^, and are scheme dependent. With the potentials used here the matching onto 



Ci(l) is known in the MS scheme at two-loop order |]T9|, This matching condition is not 

(2) 

affected by our new soft operators or the potential operators Oki^ci since their Wilson 
coefficients vanish identically at the hard scale. 
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(1) (T) 

FIG. 9: Contribution to the anomalous dimension for the production current from O^/, • 



The currents in Eq. (^) receive a non-trivial NLL anomalous dimension from graphs 
starting at order which were computed in Ref. [Q. However, this anomalous dimension 



is affected by our Oki operators which were not included there. The additional contribution 
is through the graph in Fig. R| and we find the result 



d 



167r2 



V, 



k2 \ 



(54) 



where = 2 for this spin-triplet coefficient. The anomalous dimension for the spin singlet 
coefficient Cq is identical to 7ci(^^), but with = 0. The last term in Eq. (^) which is 
proportional to a1{mv) is the contribution from Fig. ^. However, it is easy to see from 
Eq. (|5l| ) that in this case the sum of the last two terms is simply equal to Vlfl^{v)/2 so that 



167r2 



(55) 



To see how this comes about note that the counterterms for the graph involving the l/|k| 
potentials and the one in Fig. are respectively 



fc. = -v. 

5ci = — 



(56) 



k2 ' 



(^) 



cflf 



However, the factor of two difference is made up for by the fact that there is an additional 
relative factor of two in determining the anomalous dimensions i/d/duVjf^ = — 2e V^*'' + . . ., 

4e[a^V^f ] + ... in dimensional regularization (to be compared with 



(«)i 

ki 



and iyd/diy[a1Vi 

fid/dfias = —2eas + ...). This difference stems from the factors /Xg and in Eqs. (|]) 
and (|4l|) , respectively. In the next section we will consider an example where the Vki,k2 
coefficients do not come in the linear combination in Vfc,efr- 

Solving Eq. ( |5lD we find [z = as{mh')/as{m), w = as{mh''^) / as{mh')] 



In 



Ci(l)J 



-1-02 7Tas{m) {1 — z) + 03 7ras(m) ln(z) 



z 



1-13Ca/(6/3o) 



+a4 7ras[m 
+ao7ias{m) (z — 1) — w"^ ln(w) 
where the coefficients 04 and 05 agree with Ref. 



+ as Tras{m) 



z 



1-2Ca/Po 



04 



13Ca(6/5o-13Ca)2 



as 



C|[Ca(15-14S^)+/3o(4S^-3 
6(/3o-2C^)2 



(57) 



(58) 
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where = 2. The modifications to the running potentials appearing in Eq. (0) cause the 
remaining terms to differ 



02 



Cf[Ca Cf{9Ca - WOCf) - Po {^'oC\ + IQC^C^ - 32C2 )] 



as ^ ^g(6/^^_,3g)(^^_,^^) {cli^CA - IOOCf) + P. + C.(13S^ - 42) 

-6/52[2C^ + Ca(S2-3)] }, 

8C^(C^ + Cp)(C^ + 2Cp) 
ao = . (59) 

For the solution for co(z/) one should substitute = in the coefficients. 

It is important to note that the anomalous dimension 7ci('^) in Eqs. ( [5^ . |55| ) arises from 
divergences in potential loop diagrams which must be computed with a dynamic fermion 
propagator i/[E — p^/(2m)] \T^\. In Ref. [T^ it was shown that this anomalous dimension 
correctly reproduces the a^ln ag terms in the wavefunction at the origin first computed 
in Ref. For this to be the case it was necessary to include both soft and ultrasoft 

contributions to the running potentials on the right hand side of Eq. ([5^). As the scale for 
this anomalous dimension travels from m to mv it was necessary to simultaneously have 
the soft loop contributions vary from m to mv at the same rate, and have the ultrasoft 
loops vary from m to mv"^. Furthermore, in the Appendix of Ref. ||2^ it was found that 
the relation nu = n'^/m is required for a consistent subtraction of sub divergences coming 
from ultrasoft and potential divergences in three-loop vertex diagrams that contribute to the 
NNLL anomalous dimension of the production currents. This shows that ultrasoft gluons 
are needed starting at the scale m. It also shows that it is necessary to simultaneously have 
divergent loops with ultrasoft, soft and potential momentum in the theory. 

We believe that for dynamic quarks these facts are difficult to reconcile from the point of 
view of the mNRQCD-pNRQCD formalism. The existence of mNRQCD-pNRQCD seems 
to depend crucially on there being a non-trivial stage of matching that occurs at a scale 
H ~ mv ~ 1/r where soft gluons are integrated out, while the above results seem to indicate 
that such an intermediate matching scale does not exist in general. This does not mean that 
results obtained with the mNRQCD-pNRQCD formalism are necessarily wrong for dynamic 
quarks, but seems to imply that they may require some reinterpretation. 

In Ref. a procedure was developed which reproduces the anomalous dimension in 



Eq. (^) in the framework of the mNRQCD-pNRQCD formalism. In our opinion, the 
procedure suggested in Ref. contradicts some features of mNRQCD-pNRQCD upon 



which other results seem to rely. In particular, in Ref. |^ it was proposed to a) demand 
a correlation between the energy and momentum scales but only for the Ci computation, 
b) transport the 1/r matching scale back up to m by hand so that ultrasoft gluons exist 
for all scales fi < m, and c) allow couplings associated with soft loops to become unfrozen 



and run down again from m. With this construction the potentials used in Ref. [25| on the 



RHS of Eq. (|55D agree with ours, so the solution in Eq. (Q) agrees with the one in Ref. |25 



^ This result is not changed by our results here since the coefficient of the first log in the resummed V^^^ 
coefficients is the same as in Ref. and / diy/iylniy = (In^ i')/2. 
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setting Vp = mv. For this to be the case it was crucial that it was the combination V^''^^ that 



appeared in Eq. (|55D, and that the corresponding contributions in Ref. |^ were treated in 
the same way as they are predicted to be treated by vNRQCD. We agree with Ref. ||2^ that 
static quarks can be used to simphfy certain calculations, in particular those with soft loops. 
However, our conclusion is then that mNRQCD does not exist by itself as a physical theory 
that can be used to make predictions with dynamic quarks. 

Our results for the running of Ci{v) and the V^^\v) have implications for the case where 
Aqcd becomes comparable to mv"^. Take mv'^ ~ Aqcd and consider decreasing u from 
z/ = 1. As the cutoff on momentum transfers and quark momenta, mv, gets close to mw, 
the scale in the ultrasoft as{iiu) couplings is approaching the scale Aqcd- Since ultrasoft 
gluons renormalize the potentials these effects are tied together as soon as ultrasoft gluons 
first start to renormalize operators. Thus, as mv gets close to mf , a scale affecting the 
coefficients Vj(z/) of the potentials approachs Aqcd-^ Considering one ultrasoft gluon the 
affected potentials include the spin- independent l/rn? potentials, and effective l/|k| and 
potentials (or more properly the Od and O^i operators). Furthermore, due to the 
correlation between E and p^/(2m) for the heavy quarks, non-perturbative effects at the 
energy scale could lead to non-perturbative effects for the momenta. One conclusion is that 
for mv"^ ~ Aqcd we can become sensitive to non-perturbative scales through the potentials 
even for cutoffs near the momentum transfer mv ^ Aqcd- This seems quite problematic for 
perturbatively matching onto the potentials at a scale /i = mv for dynamic quarks. However, 
this does not affect the matching for static quarks as done with the mNRQCD-pNRQCD 
formalism in Ref. |]T3[. For hb states with dynamic quarks and mv'^ ~ Aqcd, this might also 
imply that non-perturbative effects have a larger influence than one usually infers. 



VIII. TOP PRODUCTION AT THRESHOLD AND QUARKONIUM ENERGIES 

The results for the order w"^, v*^ and QCD potentials and the new operators Oki and 
Oci presented in the previous sections affect the results for the NLL and NNLL top threshold 
e^e~ cross section that were given in Refs. EDI and the NNLL quarkonium energies given 



in Ref. [|r^]. Except for the handling of the Oki and l/|k| potentials discussed below, the 
changes all involve simply substituting in the coefficients and Vr^^ given in Eq. (^) and 
replacing v'f^ by the V^^^g given in Eq. (p2D. 

is) 



The modified results for the coefficients of the order spin-independent potentials, V; 



2 



is] 

and Vr , and the production current, ci, affect the results only trivially through the modified 
running of the coefficients. They do not lead to any change in the analytic form of the NLL 
and NNLL corrections to the current correlator Ai. For the corrections caused by the 
(Coulomb) potential and the operators Od it is sufficient to consider the effective Coulomb 
potential coefficient in Eq. (p^) because the corresponding matrix elements do not lead to 



UV divergences which affect the cross section. Thus, in the results presented in Ref. [0, |20 



one simply has to be replace vi*'' — >■ V^^^g and there is no other change. In Table. | we 
compare results for the coefficients at v = 0.15 and u = 0.3. The changes are smaller for 



^ The evolution towards the scale Aqcd due to ultrasoft gluons can still be computed perturbatively. This 
is very much like the fact that the evolution for the region Aqcd < /i < '\n B —t D decays can be 
computed perturbatively p^. One must just be careful not to evolve too close to the scale Aqcd- 
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matching (u = 1 


) 


-1.81 





old results {u = 


3) 


-1.72 


0.361 


new results {v = 


.3) 


-1.68 


0.359 


old results {u = 


15) 


-1.51 


0.616 


new results {v = 


.15) 


-1.39 


0.609 







c,{u)/c,{l) 


-0.357 


-1.800 


1.000 


-0.256 


-2.153 


1.034 


-0.238 


-2.153 


1.034 


-0.043 


-2.423 


1.046 


0.016 


-2.425 


1.044 



TABLE I: Comparison of the numerical change to the Wilson coefficients for ti {nit = 175 GeV). 
The old results show the coefficients used in Ref. [^], while the new results use the coefficients 
from Sects. ^ and VI. For the first three columns and the fourth column we use 1-loop running 
for as, while the last column uses 3- loop running. The coefficient vi*^(i^) is unchanged from the 
result in Ref. 



151 and is not shown. 



the larger value of u as expected. For u = 0.15 we see that the changes in Vc*"* and 
Ci(z/) are still very small (1.1%, 0.08% and 0.2%, respectively), while Vr'^^ changes by a more 
moderate amount (8%). For completeness Vj-^g is also shown, even though it is not just this 

combination of coefficients that appears in the cross section. Since V^^eff(^) ^ ^^^^ near 
u ~ 0.15 the relative change for this value of u is quite large. A more relevant measure is 

(s) 

the suppression of at ~ 0.15 compared to at = 1, which is observed in both results. 
The effect of these changes on the cross section are discussed later in this section. 

The new results for the potentials and the operators O^i lead to analytic changes in the 
NNLL corrections to the current correlator Ai because the corresponding matrix elements 
are UV divergent and the n-dependent contributions that arise from summing the potential 
indices in the operators Oki in dimensional regularization lead to modifications of the UV- 
finite terms. With dimensional regularization the full n-dependent (n = d — 1) expression for 
the potential that appears in the Schrodinger equation from the operators Oki (and Od) is 
given in App.0. Altogether, there are now three different types of corrections to the current 
correlator 5G^, (5(7^^, and 60^"^, that are needed to account for the corrections originating 
from the singlet l/|k| potential and the operators Oki and Ok2, respectively. Using MS we 
find 



a, f , m, u) 



3m^ 
4:71 a 





In ( 




^^iv — a 









+ 



3m2 
47r a 



11 

37 



17 , f ia 

hln2 + 7£; + ^ 1 

12 ' V 2f 



89 
"9 



m 



2tt a 





In ( 




^^iv — a 









+ 



m 

2-7ra 



16 



13 



21 , / ia 

— + ln2 + 7s + ^ 1- — 
16 V 2f 



175 



(60) 



In deriving these equations we have included the counterterm generated by renormalizing 
the Ji current at NLL order in Eq. (|56|). These counterterm graphs are sufficient to can- 
cel all sub divergences. The remaining overall divergences shown in Eq. (BOl) are canceled 
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by counterterms for the current correlator. For comparison the term 5G^ from the l/|k| 

is: 

In ( — 1 -2 + 21n2 + 7B 



potential which was given in Ref. |20 

IV — a 



SO (a, V, m, v) 



m 



+ 



m 



in a 



57r a 
a" 



V 



I a 

m [i 

2v 



3e2 



1 



In 2 



- -81n2 + - 
3 3 



ln^2 



TT 

"9 



(61) 



where this expression is also renormalized in the way described above. Since the 5G^'^ are 
not proportional to 5G^, it is in general not the linear combination of coefficients defined in 
V^^gg in Eq. (0) that appears in the cross section. The NNLL vector correlator 



Ai{v,m,v) = j d'^x e'^-^ (^0 T Ji,p(x) p,(0) 
p,p' 







(62) 



can be expressed in terms of the renormalized Greens functions as 



(63) 



+ Vl'\u) SG\a, V, m, v) + a^imv) V^^} {y) bG^^{a, v, m, v) 
-CFa^(mi^) V^?(z/) 56-^=^(0, w,m,z/) + bG^''\a,v,m,v) 



4 TT 



47r 



cefT 



V] 



(64) 



where q 



(-^i — 2m, 0), m is the top quark pole mass, 

^(5,r,kin/ 



= (v^ - 2m + iVt)lmfl'^, 

and the expressions for G'^{a,v,m,h'), and 6G^'^''^^^{a,v,m,h') are given in Ref. [^. The 
combination Im[c^(z/)^i] then appears in K", which is the normalized vector current induced 
cross section. From analyzing Eq. (163D we see that while the coefficients of some functions 

(s) 

are proportional to V^gg, the coefficients of other terms are not. For instance we have terms 



such as V^'^gg ^^(2), however the combination (34V^^'' + 21V^2'')^(^) ^Iso appears. Therefore, 
in general the effective l/|k| potential is not sufficient and the individual expressions for the 
Vfc,fci,fe2('^) contributions are required. 

Our results for the V-^'' coefficients also numerically affect the NNLL relation between 



is) 



the IS mass and the pole mass, which is needed to switch to the IS mass scheme This 
relation is obtained from the ^5*1 ground state solution of the Schrodinger equation and is 



given in Eq. (46) of Ref. |20|. The modifications caused by the order potentials are again 
trivial, and only the coefficients Vr*^ and V2*''are changed. For the NNLL energy levels the 
matrix elements of the Od and Oki operators do not cause additional UV divergences, so the 



effect of these operators can be implemented simply by replacing Vc'*^ and V^*'' in Ref. ||20 



by gg and vjflf^ respectively. Thus, the additional correction to the current correlator ^1 
caused by switching to the IS mass scheme has the same functional form as before 



6Ai(v,M 



IS 



A! 



^G\a,v,M 
dv 



IS 



(65) 
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8^ 
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FIG. 10: Updated numerical values of QtR" in the IS mass scheme. Also shown is a figure 
comparing the old NNLL results from Ref. [^] (solid curves) and our new results (dashed curves) 
for two values of v (red ly = 0.1, and black v = 0.15). Near the peak it is easy to see that the scale 
uncertainty is several times larger than the difference between the old and new results. 



where is the LL zero-distance Greens function of the Schrodinger equation given in 
Eq. (34) of Ref. [^. The formula for the heavy quarkonium spectrum for arbitrary quantum 
numbers is also affected by simply replacing coefficients in Eq. (40) of Ref. [|18| just as in the 
expression for the ^Si ground state energy described above. 

Next we turn to how our results numerically affect the vector-current-induced top thresh- 
old cross section K" at NLL and NNLL order. We will see that the change is quite small 
and well within the error estimate of Ref. Note that the axial- vector cross section i?" 
does not receive any modifications as it depends only on the LL value of Vc'^'^ at this order. 

In Fig. we have displayed our updated values of Q^R" in the IS mass scheme as a 
function of the cm. energy ^/s for M^^ = 175 GeV, as{mz) = 0.118 and F^ = 1.43 GeV, 
and taking four-loop running for with Uf = 5 active massless quark fiavors. The nu- 
merical methods used to evaluate Gc were described in Ref. We also show a fig- 
ure comparing our NNLL results to those in Ref. |^ for u = 0.1 and u = 0.15. The 
relative deviation is quite small and essentially independent of the cm. energy at NLL 
and NNLL order. Note that this is obvious for the NLL cross section, because at this 
order only the coefficient ci is affected by the new results. For i/i = 350 GeV the 
relative deviation amounts to (2.2%, 1.0%, 1.0%, 0.6%, 0.3%, 0.1%) at NLL order and to 
(1.5%, 0.6%, 0.2%, 0.01%, 0.1%, 0.1%) at NNLL order for u = (0.1, 0.125, 0.15, 0.2, 0.275, 0.4). 
This is well within the relative uncertainty of ±3% for the normalization of Q^K" at NNLL 
order which we estimated in Ref. |^ 



In Fig.|TT|a we have displayed Q^R" up to NNLL order. The curves show the LL (dotted 
glue lines), NLL (dashed green lines) and NNLL (solid red lines) cross section for u = 
0.1, 0.125, 0.2 and 0.4, where at LL order lower curves correspond to larger values of u. The 
conclusions drawn from these results are the same as in Ref. . Thus, compared to earlier 
NNLO (fixed order) results the variation of the normalization of the NNLL cross section 
with u is considerably reduced. Equally important, the sum of all NNLL corrections is about 
an order of magnitude smaller than the size of the NNLO corrections ||3^, and agrees well 
with the predictions at NLL order. This indicates that the expansion is converging. 

In Figs.|ll|b, c and d we have displayed separately, the contributions to Q^R" coming 
from G'^ (dashed lines), from the sum of SG^'^^'^"^ (dotted lines) and the sum of ^G'<^''~'kin,is 
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FIG. 11: Panel a) shows the NNLL results for Q'^R" with M^^ = 175 GeV at LL (dotted lines), 
NLL (dashed lines) and NNLL (solid lines) order. For each order four curves are plotted for u = 0.1, 
0.125, 0.2 and 0.4. Panels b)-d) show the relative scale dependence of various NNLL contributions 
to Q^R" for different cm. energies. The contributions are divided into those from (dashed 
lines), the sum off §Q^'^^'^'^ (dotted line), and the sum of JG^'^''^^'^'^^ (dot-dashed hues) while the 
solid (red) lines denote the sum of these terms, QtR" ■ Note that the plots have different scales for 
the y-axes. 



(dot-dashed lines) as a function of u for ^/s = 347, 350 and 353 GeV. In general we find for 
u > 0.15 that the z/-variation of the contributions from and from the sum of SG^'^^'^'^ 
cancel to a large extent, whereas the contributions from SG^'^'^^^'^^ are almost //-independent. 
On the other hand, for values of u < 0.15 and energies around and above the peak position 
the i/-variation is dominated by the contributions from SG'''^^'^'^, which rapidly increase 
in size for decreasing u, whereas the contributions from G^ and SG^'"^'^^"^'^^ are small. The 



behavior of these results are in agreement with our previous results in Refs. ||T9|, |20[ and shows 
that the modified results for the potentials and the effects of the new operators Od and Oki 
lead only to small numerical changes while the essential properties remain unchanged. To 
be more specific, making a conservative estimate by using the value 0.1 as the lower bound 
for the velocity scaling parameter to determine the remaining theoretical uncertainties 
of the NNLL renormalization group improved cross section we find 

— = ±3 % . (66) 
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This agrees with our earher estimate in Ref. and is an order of magnitude smaller than 
the uncertainties associated to fixed order NNLO QCD computations p8[. In particular, 



the conclusions that we have drawn in Ref. |2y] concerning the theoretical uncertainties in 
extractions of a^, the top Yukawa coupling and the total top width from a threshold scan at 
a future linear collider remain unchanged and are comparable to the expected experimental 
uncertainties [B^]. 



IX. CONCLUSION 



In this paper we have reconsidered the renormalization group improvement of Wilson 
coefficients in NRQCD in light of two new observations about the structure of ultrasoft 
renormalization of operators. 

We first showed that the ultrasoft renormalization of operators with soft gluons can 
induce through mixing operators O2/ whose Wilson coefficients vanish at the matching 
scale. Taking four quark matrix elements of these operators then causes a renormalization 
of the spin-independent potentials. Using the notation for the potential in 

Eq. this analysis affects the running of the coefficients V2(z^) and Vr(z^) in QCD. Our 
results are different from Ref. |15| where these additional operators were not included. We 
also compared our results to those in the mNRQCD-pNRQCD approach in Ref. p3. Our 



results do not agree with the running of 4-quark operators in mNRQCD because we find 
that the renormalization from ultrasoft gluons is present for all scales /i < m. We did find 
agreement with the pNRQCD results if we imposed a correlation between the mNRQCD- 



pNRQCD matching and renormalization scales, 



and 1/r = mu. This agreement 



is encouraging, however our correlation requirement may have implications for the NRQCD 
analysis of the case mv"^ ~ Aqcd, as discussed in Sec. |V11| and mentioned below. 

Second we performed an analysis of graphs containing mixed ultrasoft and potential loops 
and gave a new procedure for subtracting divergences in these graphs. In particular we find 
that ultrasoft gluons renormalize the operators O^i and Od displayed in Eqs. and (^4|), 
rather than the l/|k| and potentials. Because of this our results differ from those 
in Refs. |1T7| , |18|, ^ In certain situations with finite matrix elements it is possible to 
make use of effective l/|k| and 1/k^ potentials with modified Wilson coefficients. However 
in general this is not the case. An example of the former are predictions for the NNLL 
perturbative quarkonium energy levels, while an example of the latter are current correlators 

ti cross sections as discussed in Sec. 



VIII 



which give the e^e" 

We also considered the implications of the modified running for the evolution of the 
production current (Sec. [VI H ) and for predictions for the ti cross section at NNLL order. 
For the cross section the change from the results in Ref. |19|, ^ is very small, being ^ 1% 
for the physically motivated values u > 0.15. If we include a larger more conservative 
range of scale variations then the change is still small becoming 1.5% at u = 0.1. Since 
the uncertainty assigned to predictions in Ref. was ±3%, all analyses and conclusions 
for the cross section are unchanged. In particular, there is still a large improvement in the 
convergence of the perturbation theory over not summing the logarithms, and we can assign 
a conservative overall uncertainty of ±3% to the NNLL cross section predictions which is 
much smaller than the uncertainty found at NNLO P^ . 

Our results for the running have implications for cases where Aqcd becomes comparable 
to mf ^. Essentially, since ultrasoft gluons renormalize spin-independent potential operators 
these potentials become sensitive to the scale Aqcd at an earlier stage than might otherwise 
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be thought. The hst of affected potentials includes the spin-independent potentials 
and the Od and Oki operators (which give effects often attributed to the l/|k| and 
potentials). Due to the correlation for mv"^ ~ Aqcd we are sensitive to non-perturbative 
scales through the potentials even though the momentum transfer mv ^ Aqcd- This seems 
problematic for perturbatively matching onto the potentials at a scale /x = mv for dynamic 
quarks, because at this scale the coefficients of the potentials are already blowing up. 
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APPENDIX A: SUMMARY OF STRUCTURES FOR SIX FIELD OPERATORS 

The structures appearing in the QCD operators 0^2!p^ ^2A^ ^^'^ ^2c include 



r(0) 



A,fj.u 



r(o) 

Alp 

p(0),CD 

c 

r(0) 



yAyB (2g°7° + k ■ 7) ^^^^ 



1 



k2 + 2k ■ q 

T-(o) ^ 



(2^070 - k ■ 7) 
k2 - 2k ■ q 



7(0) 
^0 



A_EC rBED 



1 



AED /SEC 



t^°ng,P,p') (t/(°))/(-g,-p,-p') 

t/Sn-g,p,p') (f/^°^)/(?, -p,-p') 



.(0) 

A,x 



jAEC jBED jAED jBEC 

k2 + 2k ■ q ^ k2 - 2k ■ q 

A p(0) _ 'fB 



(0) 



1 



k2 + 2k ■ q 

1 

k2 - 2k ■ q 



(Al) 



where the U^S) ^ \ ^-nd Y^^^ coefficients can be found in Ref. [1^. We have also made use 
of color contractions of these structures including 



p(0),(l) ^ p(0),(T) ^ ^ 
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^{0),(T) 
A,/ii^ 
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.(o),(r) 



r(o),(T) 



(2^07° + k ■ 7) ^ (2907° - k ■ 7) 
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APPENDIX B: FORM OF POTENTIALS IN THE SCHRODINGER EQUATION 



When the l/|k| potential and the operators Oki and 0^2 are included in the Schrodinger 
equation in dimensional regularization the form of the potential is 



ySchr. 



p.q) 



-7r2 ,,2e 



2^2- 



m(k' 



^ alimu) V^IV) - Sri - 12) /(1, 1) 



k2 ^ 



— CpaJmi') V, 



4(n2-3n + 8)/(l,l) 



m(k2)^" 2 

while for the 1/k^ potential and the operators Od, Oc2, and Cc3 the potential is 

^3,,6e 



(Bl) 



Schr. / 



lP,q) = ^^i'\^) + ^^CFasimvf 



n 



where 



/(a, 6) 



(4-n)/(l,2)/(3--,l) 



r(a + &-f)r(f -a)r(f -6) 

r(a)r(6)r(n-a-6)(47r)t 



(B2) 



(B3) 



In Eq. ( P2| ) the operator Od gives a vanishing contribution. Note that the functions 
l/(k^)^~'^^ and l/(k^)^~" are not the d-dimensional Fourier transform of a or 1/r 
potential. Also note that in another regularization scheme such as with a cutoff, the poten- 
tials in Eqs. ( pi|) and ( P2|) would take a different functional form, however the momentum 
dependence of the Vjt, Vki,k2, Vc, and Vc2,c3 terms would still differ from each other. 
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